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1. Introduction 



One of the classical problems of mathematics consists in finding a closed 
plane curve of a given length that encloses the greatest area: the isoperimetric 
^ ' problem. The legend says that the first person who solved the isoperimetric 

problem was Dido, the Queen of Carthage, who was offered as much land as she 
could surround with the skin of a bull. Dido's problem is nowadays part of the 
calculus of variations [23j[35] . 

Fractional calculus is a generalization of (integer) differential calculus, 
allowing to define derivatives (and integrals) of real or complex order |25U301l32] . 
The first application of fractional calculus belongs to Niels Henrik Abel (1802- 
1829) and goes back to 1823 [lj. Abel applied the fractional calculus to the 
solution of an integral equation which arises in the formulation of the tautochrone 
problem. This problem, sometimes also called the isochrone problem, is that of 
finding the shape of a frictionless wire lying in a vertical plane such that the 
time of a bead placed on the wire slides to the lowest point of the wire in the 
same time regardless of where the bead is placed. The cycloid is the isochrone 
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as well as the brachistochrone curve: it gives the shortest time of slide and marks 
the born of the calculus of variations. 

The study of fractional problems of the calculus of variations and respec- 
tive Euler-Lagrange type equations is a subject of current strong research due 
to its many applications in science and engineering, including mechanics, chem- 
istry, biology, economics, and control theory [27]. In 1996-1997 Riewe obtained 
a version of the Euler-Lagrange equations for fractional variational problems 
combining the conservative and nonconservative cases [33,34J. Since then, nu- 
merous works on the fractional calculus of variations, fractional optimal control 
and its applications have been written— see, e.g., [Ill7llMTT HMI201^l26ll2B] and 
references therein. For the study of fractional isoperimetric problems, see [5]. 

In the pioneering paper [2], and others that followed, the fractional nec- 
essary optimality conditions are proved under the hypothesis that admissible 
functions y have continuous left and right fractional derivatives on the closed 
interval [a, b] . By considering that the admissible functions y have continu- 
ous left fractional derivatives on the whole interval, then necessarily y(a) = 0; 
by considering that the admissible functions y have continuous right fractional 
derivatives, then necessarily y(b) = 0. This fact has been independently re- 
marked, in different contexts, at least in [5|ITT | [T3 | [2l] . 

In our work we want to be able to consider arbitrarily given boundary 
conditions y(a) = y a and y(b) = y\, (and isoperimetric constraints). For that we 
consider variational functionals with integrands involving not only a fractional 
derivative of order a G (0, 1) of the unknown function y, but also the classical 
derivative y'. More precisely, we consider dependence of the integrands on the 
independent variable t, unknown function y, and y' + k a Dfy with k a real pa- 
rameter. As a consequence, one gets a proper extension of the classical calculus 
of variations, in the sense that the classical theory is recovered with the partic- 
ular situation k = 0. We remark that this is not the case with all the previous 
literature on the fractional variational calculus, where the classical theory is not 
included as a particular case and only as a limit, when a — > 1. 

The text is organized as follows. In Section 2 we briefly recall the neces- 
sary definitions and properties of the fractional calculus in the sense of Riemann- 
Liouville. Our results are stated, proved, and illustrated through an example, 
in Section 3. We end with Section 4 of conclusion. 

2. Preliminaries 

In this section some basic definitions and properties of fractional calculus 
are given. For more on the subject we refer the reader to the books 
and historical survey [27] , 
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Definition 1. (Left and right Riemann-Liouville derivatives) Let / be 
a function defined on [a, b}. The operator a Df , 

a D*f(t) = - 1 D n [\t - rr-^ l f(r)dr , 
T{n - a) J a 

is called the left Riemann-Liouville fractional derivative of order a, and the 
operator jD", 

T{n - a) J t 

is called the right Riemann-Liouville fractional derivative of order a, where 
a E M + is the order of the derivatives and the integer number n is such that 
n — 1 < a < n. 

Definition 2. (Mittag-Leffler function) Let a,/3 > 0. The Mittag- 
Leffier function is defined by 

Theorem 3. (Integration by parts) If f,g and the fractional derivatives 
a Dfg and tD^f are continuous at every point t S [a, b], then 

f f(t) a D?g(t)dt = f git) t Dlfit)dt (1) 

J a J a 

for any < a < 1. 

Remark 4. If f(a) / 0, then a Dtf(t)\ t=a = oo. Similarly, if fib) / 0, 
then tD^f(t)\ t=b = oo. Thus, if / possesses continuous left and right Riemann- 
Liouville fractional derivatives at every point t £ [a,b], then f{a) = fib) = 0. 
This explains why the usual term f(t)g(t)\ a does not appear on the right-hand 
side of CD). 

3. Main results 



Following [24J, we prove optimality conditions of Euler-Lagrange type 
for variational problems containing classical and fractional derivatives simulta- 
neously. In Section 3.1 the fundamental variational problem is considered, while 
in Section 3.2 we study the isoperimetric problem. Our results cover fractional 
variational problems subject to arbitrarily given boundary conditions. This is 
in contrast with [2 H4"lll5| , where the necessary optimality conditions are valid for 
appropriate zero valued boundary conditions (cf. Remark 4). For a discussion 
on this matter see 
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3.1. The Euler Lagrange equation 



Let < a < 1. Consider the following problem: find a function y £ 
C 1 [a, b] for which the functional 

J(y) = f F(t, y(t),y'(t) + k a D?y(t)) dt (2) 

J a 

subject to given boundary conditions 

y(a) = Va, y{b) = yb, (3) 

has an extremum. We assume A; is a fixed real number, F G C 2 ([a,b] x R 2 ;R), 
and d^F (the partial derivative of F(-,-,-) with respect to its third argument) 
has a continuous right Riemann-Liouville fractional derivative of order a. 

Definition 5. A function y £ C l [a, b] that satisfies the given boundary 
conditions (J3J) is said to be admissible for problem (|2"])-( f 5 ]) . 

For simplicity of notation we introduce the operator [•]£ defined by 
[y] a k (t) = {t,y(t),y'(t) + k a D?y(t)) . 
With this notation we can write ([2]) simply as 

J(y)= [ F[y]t(t)dt. 

J a 

Theorem 6. (The fractional Euler-Lagrange equation) If y is an 
extremizer (minimizer or maximizer) of problem ([2|)-([3]) ; then y satisfies the 
Euler-Lagrange equation 

d 2 F[y]Ut) - j t d,F[y]t{t) + k t D^d 3 F[y] a k (t) = (4) 

for all t £ [a, b] . 

Proof. Suppose that y is a solution of ©-Q. Note that admissible 
functions y can be written in the form y(t) = y(t) + erj(t), where rj £ C 1 ^,^], 
77(a) = 77(6) = 0, and e £ R. Let 

J(e) =J a F U V(t) + «?(*). J f (V(t) + £V(t)) + W (y(t) + erj(t))\ dt. 

Since a -Df is a linear operator, we know that 

a D« (y(t) + e V (t)) = a D?y(t) + e a D? V (t). 
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On the other hand, 



dJ 
~dl 



e=0 J a 



j/[y\ a k {t)dt 



<E = 



drj(t) 



d 2 F[y] a k (t) ■ V (t) + d3F[y]%(t)-^- + kd 3 F[y}<Z(t) a D?r)(t))dt. 



Using integration by parts we get 



and 



/«Ff*-*^i-j((4*F)* 



f 8 3 F a D?ridt= f r,,Dld 3 Fdt. 

a J a 



(5) 



(6) 



(7) 



Substituting ([6]) and $7$) into (jSJ), and having in mind that i](a) = rj(b) = 0, it 
follows that 



dJ_ 
~d~e 



e=0 J a 



d 



v(t)(d 2 F[y] a k (t) - -d 3 F[y] a k (t) + k t D a b d 3 F[y]t(t))dt. 



A necessary optimality condition is given by ^fi _ = 0. Hence, 



ri{t)[d 2 F[y\t(t) - -d 3 F[yr k (t)+k t D^d 3 F[y}^t))dt = 0. (8) 



We obtain equality §%$) by applying the fundamental lemma of the calculus of 
variations to (jSJ). ■ 

Example 7. Note that for k = our necessary optimality condition 
dH) reduces to the classical Euler-Lagrange equation 



3.2. The fractional isoperimetric problem 

As before, let < a < 1. We now consider the problem of extremizing a 
functional 

J(y) = I F(t, y(t),y'(t) + k a D?y(t)) dt (9) 

J a 

in the class y £ C 1 ^, b] when subject to given boundary conditions 

y(a)=Va, y{b) = y b , (10) 
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and an isoperimetric constraint 

X{y) = f G(t,y(t),y'(t) + k a D?y(t))dt = Z. (11) 

J a 

We assume that k and £ are fixed real numbers, F,G £ C 2 ([a,b] x M 2 ;R), and 
d 3 F and c^G have continuous right Riemann-Liouville fractional derivatives of 
order a. 

Definition 8. A function y £ C 1 [a, 6] that satisfies the given boundary 
conditions (jlOp and isoperimetric constraint (jlip is said to be admissible for 
problem pi-dTTT). 



Definition 9. An admissible function y is an extremal for X if it satisfies 
the fractional Euler-Lagrange equation 

d 2 G[y} a k (t) - j t d 3 G[y\Ut) + k t D^d 3 G[y] a k (t) = 

for all t £ [a, b]. 

The next theorem gives a necessary optimality condition for the fractional 
isoperimetric problem (j9l)~ (|lip . 

Theorem 10. Let y be an extremizer to the functional ([9]) subject to 
the boundary conditions (jlOp and the isoperimetric constraint (jlip . If y is not 
an extremal for X, then there exists a constant A such that 

d 2 H[y]Ut) - j t d 3 H[y}%(t) + k t D%d 3 H[y}Ut) = (12) 

for all t £ [a, b] , where H(t, y , v) = F(t, y , v) — XG(t, y, v) . 
Proof. We introduce the two parameter family 

V = V + ei?7i + e 2 ?72, (13) 

in which rj\ and r\2 are such that 7/1,7/2 £ C x \a % b] and they have continuous left 
and right fractional derivatives. We also require that 

7/1 (a) = 7/1(6) = = 772(a) = 7/2(6). 

First we need to show that in the family (|13|) there are curves such that y satisfies 
(fTTj) . Substituting y by y in (fTT|) . X(y) becomes a function of two parameters 
ei,e 2 . Let 

rb 



/(ei,e 2 )= / G(t,y,y' + k a D?y)dt-Z- 

J a 



Calculus of Variations with Classical 
Then, 1(0, 0) = and 



de 2 



(0,0) 



J m (d 2 G - j t d 3 G + ktDSdaG) dt. 



Since y is not an extremal for X, by the fundamental lemma of the calculus of 
variations there is a function 7] 2 such that 



ai_ 

de 2 



^0. 



(0,0) 



By the implicit function theorem, there exists a function e 2 (-) defined in a neigh- 
borhood of zero, such that I{e\, £2(61)) = 0. Let J(e±, e 2 ) = J{y)- Then, by the 
Lagrange multiplier rule, there exists a real A such that 

V(J(0,0) -A/(0,0)) = 0. 



Because 



dJ_ 



(0,0) 



j m (d 2 F - j t 3 F + k t D%d 3 F) dt 



and 






(0,0) 



J Vi I.O2G - j f d 3 G + k t D%d 3 G) dt, 



one has 



m 



d 2 F - j t 8 3 F + k t D%d 3 F) - A (d 2 G - j t d 3 G + k t D%d 3 G 



dt = 0. 



Since r\\ is an arbitrary function, (|12|) follows from the fundamental lemma of 
the calculus of variations. ■ 
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3.3. An example 

Let a G (0, 1) and k, £ G R. Consider the following fractional isoperimet- 
ric problem: 

.76/)=/ (y'+^o^) 2 A^min 
JO 

z(y) = / {l/ + k I%y)dt = Z (14) 

JO 

y(0) = , 2/(1) = y ^ (-& (1 - r) 1 " ) £dr. 

In this case the augmented Lagrangian H of Theorem 10 is given by H(t, y, v) = 
v 2 — Xv. One can easily check that 

y(t) = J £ w (-* (t - r) 1 -") ^dr (15) 

• is not an extremal for X; 

• satisfies y' + k^Dfy = £ (see, e.g., [251 P- 297, Theorem 5.5]). 
Moreover, (JT5J) satisfies CE]) for A = 2£, i.e., 

-^ (2 (i/ + k D?y) - 2^ + k t D? (2 (y' + k D?y) - 2£) = 0. 



We conclude that (|15|) is the extremal for problem ([T] 

Example 11. Choose k = 0. In this case the isoperimetric constraint 
is trivially satisfied, (|14p is reduced to the classical problem of the calculus of 
variations 



J(y) = / (y'(t)) 2 dt 



mm 



/o (16) 

y (0) = 0, y(l)=£, 

and our general extremal f)15j) simplifies to the well-known minimizer y(£) = £t 
of (USD. 



Example 12. When a — >• 1 the isoperimetric constraint is redundant 
with the boundary conditions, and the fractional problem (|14p simplifies to the 
classical variational problem 

J(y) = {k + l) 2 [ y\t) 2 dt — ► min 

^ (17) 

y(0) =0, 2/(1) = ^-. 
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Our fractional extremal ([T5]) gives y(t) = -A^t, which is exactly the minimizer 
of (HZD- 



Example 13. Choose k = £ = 1. When a — > one gets from (|14p the 
classical isoperimetric problem 



J(y)= f\y'(t)+y(t)) 2 dt 
Jo 



mm 



;is) 



f 1 i 

X y = / y(t)dt = - 

1/(0) =0, y(l) = l--. 

e 

Extremal (I15p is then reduced to the classical extremal y(t) = 1 — e - ' of f\l8\\ . 

Example 14. Choose k = 1 and a = ^- Then (J14D gives the following 
fractional isoperimetric problem: 

J(y)=f (y' + oDfyj dt^ mm 

I (y) = J Q (y + oD?y\dt = t (19) 



2/(0) = 0, |/(1) = -e ( 1 - erfc(l) + 2 



where erfc is the complementary error function. The extremal (|15p for the 
particular fractional problem (1191) is 

2Vt\ 



l(t) =-Hl- e'erfc(v / t) + 



4. Conclusion 

Fractional variational calculus provides a very useful framework to deal 
with nonlocal dynamics in Mechanics and Physics [6l[T6]. It has received con- 
siderable interest in recent years, with several researchers applying this field to 
develop fractional classical and quantum mechanics [101 118 ( 129]. Motivated by 
the results and insights of [5,8,24], in this paper we generalize previous fractional 
Euler-Lagrange equations by proving optimality conditions for fractional prob- 
lems of the calculus of variations where the highest derivative in the Lagrangian 
is of integer order. This approach avoids the difficulties with the given boundary 
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conditions when in presence of the Riemann-Liouville derivatives [24J . For the 
case with the Caputo fractional derivatives ( [32]) we refer the reader to [31J. 

We focus our attention to problems subject to integral constraints (frac- 
tional isoperimetric problems), which have recently found a broad class of im- 
portant applications [9lll7jll9| . For k = our results are reduced to the classical 
ones [35] . This is in contrast with the standard approach to fractional variational 
calculus, where the integer-order case is obtained only in the limit. 
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